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Abstract. We have studied the low speed fracture regime for different glassy 
materials with variable but controlled length scales of heterogeneity in a carefully 
mastered surrounding atmosphere. By using optical and atomic force microscopy 
(AFM) techniques we tracked in real-time the crack tip propagation at the nanometer 
scale on a wide velocity range (10~^ — 10~^^ m/s and below). The influence of the 
heterogeneities on this velocity is presented and discussed. Our experiments reveal 
also - for the first time - that the crack progresses through nucleation, growth and 
coalescence of nanometric damage cavities within the amorphous phase. This may 
explain the large fluctuations observed in the crack tip velocities for the smallest values. 
This behaviour is very similar to what is involved, at the micrometric scale, in ductile 
fracture. The only difference is very likely due to the related length scales (nanometric 
instead of micrometric). Consequences of such a nano-ductile fracture mode observed 
at a temperature far below the glass transition temperature, Tg, in glass is finally 
discussed. 
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1. Introduction 

Despite many works and recent progresses, materials fracture still exhibits many 
puzzling aspects. Two different observations seem particularly difficult to reconcile: 

• Brittle materials - whose most common example is glass - break abruptly, without 
first deforming, while in ductile materials - like metallic alloys - fracture is preceded 
by large plastic deformations. In this latter case, the crack is observed |Ij to 
progress through the coalescence of micrometric damage cavities nucleated from 
microstructural defects (second phase precipitates, grain boundaries...). 

• However, quantitative studies reveal that both types of fracture surfaces have a very 
similar morphology. Fracture surfaces have been shown to be self-affine objects 0E] 
for both brittle and ductile materials. Two self affine regimes coexist : at small 
length scales (up to a length (^c), the roughness exponent C is close to 0.5, whereas 
at large length scale (from to a larger length ^), C, is around 0.8. These values 
of C are universal, but as well as ^ depend on the kind of material considered, 
ranging from a few nanometres for glass |lj to a few centimetres for concrete and 
rocks [5 . It is worthy to note that, for ductile materials, is of the order of the 
typical size of the damage cavities when they coalesce PJ. 

This strongly suggests that the observed transition between the two self-affine 
regimes reflects a change from an intra-cavity structure to an inter-cavity one. Such 
a scenario was first proposed in [110111, based on the observation of the roughness 
of growing cavities before coalescence in an aluminium alloy and further developed 
in a recent paper [Hj. 

How to reconciliate these two observations? A natural supposition is that a ductile- 
like - up to now undetected - fracture process exists in brittle materials like glass. In 
this case, the size of the associated damage cavities at coalescence is expected to be 
equal to the crossover length ^c, ^-e. of the order of ten nanometers. The need for an 
experimental observation of such effect was at the genesis of the experiment described in 
the following. Our new experimental set-up, based on Atomic Force Microscopy (AFM), 
is described in section 2. Huge velocity fluctuations in the progression of the tip were 
detected (section 3) and explained by the observation of nanometric cavities developing 
ahead of the crack tip (section 4). These cavities were proved to be damage cavities 
by using two independent methods (section 5). Finally, the possible origins of these 
cavities are conjectured and the consequences of such ductile fracture mode in glass are 
discussed (section 6). 

2. Experimental set-up 

The design of our experiments should meet several successive requirements to achieve 
the goal described above. Those are listed below: 
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Figure 1. Experimental setup. 

(i) The apparatus should have a very high spatial resolution (a few tenths of 
nanometers) and work in real space and real time (no post-mortem studies). 
Moreover, it should deal with electrically non conductive materials. Surface 
treatments - which limit the spatial resolution (by metal deposition for instance) - 
are also prohibited. The AFM technology appears then to be the ideal tool that 
meets all these requirements. Due to the nanometric resolution of AFM we measure 
crack speeds lower than 10~^^ m.s~^ with a reasonable acquisition time (15 minutes 
between each acquired image). 

(ii) The main drawback of AFM resides in its recording time, that cannot be set lower 
than a few minutes without significant loss in the spatial resolution. It forces us 
to work in the so-called 'stress corrosion' or 'sub-critical growth' regime Pl lim lTT|. 
The crack behaviour is then controlled by the chemically activated processes with 
the surrounding environment. A rigorous control of atmosphere and temperature 
is then necessary. 

(iii) The crack progression should be stable and well monitored. This leads us to design 
a mechanical DCDC (Double Cleavage Drilled Compression) system |12j . 

The experimental set-up is illustrated in Fig. ^ All the experiments are performed 
at a constant temperature of 22.0±0.5°C in a leak-proof chamber under an atmosphere 
composed of pure nitrogen and water vapour at a relative humidity level of 42 ± 1% 
after preliminary out-gassing. 

Lithium alumino-silicate glass-ceramics were studied. The main advantage of this 
kind of materials is that its structure can be controlled through an applied thermal 
treatment ^3]- Two types of structures are investigated below. The first one (sample 
A) corresponds to a pure glassy state as identified by X-ray diffraction (Fig. \^). A 
single thermal treatment at 660° C is performed in order to remove residual stresses. 
The second one (sample B) is classified as a slightly unglassy state in which small 
crystals of /5-quartz phase have formed (Fig. Eb) in the glassy phase. The structure 
of sample B is obtained after a two-step thermal treatment at plateau temperatures of 
Ti=750°C and T2=900°C. Both the size and the volume fraction of crystalline grains 
are evaluated by imaging the samples surface by AFM after a stay in fiuorhydric acid 
(concentration ranging from 0.4% to 2%) for 30 s, which dissolves the amorphous phase 
faster than the crystalline one pHj . 
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Figure 2. X-ray diffraction patterns on lithium alumino-silicate glass-ceramics, (a): 
Sample A in a pure glassy state, b): Sample B in a slightly unglassy state. 
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Figure 3. Sketch of the DCDC geometry. 



Fractures are then performed on the DCDC set-up |T2l US]: Parallelepipedic 
(4x4x40 mm^) samples (Fig.E)) are designed with a cyhndrical hole (radius a = 0.5 mm) 
drilled in the centre and perpendicularly to the 4 x 40 mm^ surface (Fig. IH)). The hole 
axis defines the z-direction. The x-axis (resp. y-axis) is parallel to the 40 mm (resp. 
4 mm) side of the 4 x 40 mm^ surfaces. 

In all cases, the 4 x 40 mm^ surfaces are optically polished (the measured RMS 
roughness is around 0.25 nm for a 10 x 10 fim^ scan size). A compressive load is applied 
perpendicularly to the 4x4 mm^ surfaces. The external stress a is gradually increased by 
the slow constant displacement (0.02 mm/min) of the jaws of the compressive machine 
(Fig. Q). Once the two cracks are initiated - symmetrically to the hole axis -, the jaws 
displacement is stopped. The crack then propagates along the x-axis in the symmetry 
plane of the sample parallel to the {x,z) plane. In this geometry, the stress intensity 
factor Kj is given by [T3]: Kj = a^faj (0.375c/a-|-2), where c is the crack length (Fig.Ej). 
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The crack motion within the {x,y) surface is monitored by our experimental system 
combining optical microscopy and AFM. Optical image processing gives the position of 
the crack tip (with a spatial resolution of 5 micrometers) and consequently the velocity 
for V ranging from 10~^ to 10~^ m.s~^. Let us note that the measured "instantaneous" 
velocity is actually averaged over a temporal window the width of which is given by time 
taken by the crack tip to advance on a distance equal to the spatial resolution (typically 
5000s for a mean speed of 10~^ m.s"^). By AFM measurements ^H] - performed in a high 
amplitude resonant mode ("tapping" mode) -, one probes the crack tip neighbourhood 
at magnifications ranging from 75 x 75 nm^ to 5 x 5 /xm^. The crack tip can be tracked 
at velocities ranging from 10"^ to 10^^^ m.s^^. In this latter case, the minimum width 
of the temporal window is reduced by two orders of magnitude thanks to the increased 
spatial resolution. 

3. Crack tip velocity 

At the very first moments, the crack propagates very quickly. In this regime, the crack 
velocity v is independent of the chemical composition of the surrounding environment [H]. 
As the crack length c increases, Ki decreases, and v decreases quickly. Under vacuum, 
the crack stops for Kj smaller than a critical value Ki^ referred to as the toughness of 
the material. But in a humid atmosphere, the corrosive action of water on glass allows 
for slow crack propagation at much lower values of the stress [H]. When Kj becomes 
smaller than the fracture toughness Ki^ (i.e. in the stress corrosion regime), the crack 
motion is slow enough to be monitored by our experimental system. 

Figure |^ shows the variation of the velocity f as a function of the stress intensity 
factor Ki for both the amorphous and the partially devitrified specimens. The 
observed exponential behaviour is compatible with a stress enhanced chemical activated 
process. CHI ■ 

The velocities measured for the partially devitrified sample (Fig. EJd) are shifted 
toward positive Kj when compared to the data relative to the amorphous sample. 
To understand this shift, we probe ^1] the surface crack path at sub- micro metric 
scales (Fig. EI). In the partially devitrified specimen, the crystalline germs deflect 
the crack (Fig. Eb). As the crack tip keeps going through the amorphous phase, the 
stress enhanced chemical activated process is identical in both samples (A and B) and 
consequently the slope of the semilogarithmic curves v{Kj) are the same in both the 
amorphous and the partially devitrified specimens. However, the deflections of the crack 
by the crystals induce local mode II and mode III components in the local stress intensity 
factor, which toughens the material Consequently, Kjc is larger in the devitrified 
sample, which shifts the v{Ki) curve. A similar effect of the size of heterogeneities has 
been observed for the fatigue of metallic alloys [TH]. 

The lowest velocities - only reachable thanks to very high spatial resolution of 
AFM- shown in Fig. |3] exhibit important fluctuations - of the order of the average 
velocity - for 6ot/i specimens jlH]. Consequently, they cannot be related to the crystalline 
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Figure 4. Variation of the crack tip velocity v versus the stress intensity factor Kj. 
Open circles (resp. black triangles) correspond to optical measurements (resp. AFM 
measurements). The error bars correspond to the standard deviation on the velocity 
for a fixed value of the stress intensity factor Kj. For the lowest velocities the 
fluctuations are of the order of the average velocity. 
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Figure 5. Crack path in (a) sample A (the amorphous specimen) and (b) sample 
B (the partially devitrificd one). Frames are phase images sensitive to local 
mechanical properties, which allows to distinguish the crystalline heterogeneities. In 
the partially devitrified sample, the crystals deflect the crack which propagates within 
the amorphous phase. 



heterogeneities but are more likely inherent to the amorphous phase. The velocity 
fluctuations shown in figure 4 are of the order of 5.10"^^ m/s, i.e. of the same order 
of the magnitude as the lowest measured averaged velocities : this fiuctuations are the 
same in both specimens. 

4. Evidence of nanometric cavities ahead of the crack tip 

To understand the origin of the velocity fiuctuation, we probed the neighbourhood of 
the crack tip at the nanometer scale in the amorphous specimen (Fig. IHI). This clearly 
reveals the presence of cavities of typically 20 nm in length and 5 nm in width ahead of 
the crack tip (Fig. |3i) [201 ■ 

These cavities grow with time (Fig. IHb) until they coalesce (Fig. IHt). At these 
nanometric scales, the crack front does not propagate regularly, but intermittently 
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Figure 6. Sequence of successive topographic AFM frames showing the crack 
propagation at the surface of the specimen. The scan-size is 75 x 75 nm^ and the 
heights range over 2 nm. The recording time for one frame is around 3 mn and two 
successive frames are separated by 20 mn. The crack front propagates from the left to 
the right (x-direction) with a mean velocity v of around 10~^^ m/s. (a): evidence of 
nanometric damage cavities before the fracture advance, (b): growth of the cavities, 
(c) : the crack is advancing via the coalescence of all the cavities 

through the merging of the nano-scale cavities, which explains the large fluctuations 
observed for the lowest velocitites (Fig. HJ. 

5. Nanometric cavities and nano-ductility. 

To ensure that the spots observed ahead of the crack tip are actually damage cavities 
which grow further and coalesce with the main crack leading to failure, we use the 
Fracture Surface Topography Analysis (FRASTA) technique flrst introduced to study 
damage in metallic alloys |2I1I22]- These authors showed that, in a ductile scenario, 
each cavity initiation is accompanied by local irreversible plastic deformations printed 
in relief on the developing fracture surfaces (the crack lines when the method is applied 
in two dimensions as in the present case) that should remain visible after the cavities 
have coalesced and the crack has crossed. The FRASTA method is designed to analyse 
the mismatch between the lower fracture surface (the lower crack line when the method 
is applied in two dimensions as in the present case) and the upper one (the upper crack 
line), mismatch due to the presence of cavities. For that purpose, the crack lines are 
flrst determined by binarising the image of the sample after fracture (Fig. [2K), and the 
unbroken material is reconstituted virtually by placing numerically the lower line over 
the upper one (Fig. [3d)- By translating the lower crack line gradually in the direction of 
decreasing y (Fig.[7|D), - this displacement is what actually occurs during fracture -, one 
can see the cavities appearing and growing in the chronological order. The structure 
obtained for a given displacement, i.e. at a given time, is superimposed on images 
recorded prior to failure and shown to correspond actually to cavities observed at this 
given time (Fig. [Zb)- This clearly indicates that the spots shape is determined well 
ahead of the crack tip, before the flrst ones can be actually observed, which provides a 
rather strong argument to relate these nano-scale spots to damage cavities. 
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Figure 7. Fracture Surface Topographical Analysis (FRASTA). (a) Frame 0); (broken 
sample) is binarised and the contours of the crack are determined, (b) The lower 
line is first numerically raised over the upper one and then gradually displaced in the 
direction of decreasing y, as schematised by the arrow. Cavities are coloured in green, 
(c): Result of the method: superimposition of the obtained cavities on the image Els 
recorded prior to complete failure. 



This nano-scale ductility is also confirmed by the study of the displacement field 
around the crack tip (see also references |23 121] for related discussion): for a slit-like 
plane crack in an ideal Hookean continuum solid, all the components of the stress tensor 
are predicted to decrease as r^^/^ with the distance r from the crack tip [2S]- If the 
macroscopic stress-strain relation work on the nanometer scale the depth 5 (see the 
caption of Fig. |H1 for a precise definition) should thus decrease as r~^/^ in the vicinity 
of the crack tip during crack propagation. To be more precise, 5 = Ar~^^'^ cos{9/2) for 
mode I fractures [2^] where A is a pre-factor depending only on the applied load and on 
the specimen geometry, and 6 is the angle between the direction of crack propagation and 
the r direction. Measurements of depth profiles have been performed on 1 x 1 /im^ AFM 
topographical frames (Fig. |HK) along the direction of crack propagation (Fig. |Ht>) and 
perpendicularly to it (Fig. |Ht). For both profiles, 5 departs from the linear elastic r~^/^ 
scaling for r smaller than a threshold Vc highly dependent on 6: for 6 = 0°, Vc = 100 nm 
while for 6 = 90°, = 20 nm. These short range departures from the linear elastic 
behaviour may be related to the presence of cavities although other phenomena could 
be responsible for this discrepancy (as for instance viscous effect, chemical effect, ...). 
However, the fact that the order of magnitude of the ratio rc{6 = 0°)/rc{9 = 90°) - much 
higher than the ratio of the cosine terms in the linear elastic expression of 5 - is close to 
the aspect ratio of the observed damage cavities, strongly suggests a correlation between 
damage and non linear elasticity. 

6. Discussion 

Same fracture experiments performed in amorphous silica specimens reveal similar 
damage cavities. This suggests that their existence does not depend on the precise 
chemical composition of the studied glass. The origin of the nucleation of cavities should 
be found more likely in the amorphous structure, which contains inherent atomic density 
fluctuations at the nanometer scale. Such a scenario was indeed predicted by Molecular 
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Figure 8. Measurements of the surface deformations and comparison with the 
predictions for an ideal Hookean material. The crack propagates from left to right (x 
positive), (a): typical AFM topographical frame of the vicinity of the crack tip. The 
scan-size is 1 x 1 /xm^ and the heights range over 3 nm. The white vertical (respectively 
horizontal) dotted line sets the x-coordinate Xf (respectively the y-coordinate yf) of 
the crack tip. (b) (respectively (c)) top: Plot of the z profile along (respectively 
perpendicularly to) the direction of crack propagation. The open circles correspond 
to experimental data while the full line corresponds to the prediction z = zq — Ar~^/'^ 
(where zq and A are fit parameters) given for an ideal Hookean solid. Bottom: Log-log 
plot of the depth 5 — zq — z versus the distance r = x — x f (respectively r — y — yf) 
from the crack tip. For r < rc, the (5-profile departs from the predictions given by the 
linear elastic theory. 

Dynamics Simulation [23 ESI 1201 IHO] that evidenced atomic density fluctuations in the 
structure of simulated amorphous silica : The Si and O atoms are shown to form 
silica tetrahedra connected together to build rings of different sizes ranging from 3 
to 9 tetrahedra. At larger length scales, ranging from 1.5 nm to 6 nm, the density 
of these rings is found to fluctuate with high density areas surrounded by low density 
areas. Moreover, the Molecular Dynamics simulations of van Brutzel fl^ 123 IHO] show 
that, at this length scale, crack propagates by growth and coalescence of small cavities 
which appear in areas with low density of rings, ahead of the crack tip. They behave 
as stress concentrators and grow under the stress imposed by the presence of the main 
crack to give birth to the cavities actually observed in the AFM frames. 

The comparison of these theoretical results and the experimental evidences reported 
here strongly suggests that the crack advance in glass occurs by a ductile-hke process 
- at a temperature far below Tg - bound to nucleation, growth and coalescence of damage 
cavities. The nucleation of these damage cavities is initiated very likely within the low 
density nanometric areas. The origin of the universality in the existence of two roughness 
regimes observed on post-mortem glass fracture surfaces |S1E] can now be explained: The 
fracture advances through the nucleation, growth and coalescence of damage cavities in 
both brittle and ductile material. The only difference resides in the length scales of the 
damage cavities : those related to 'brittle' materials are approximately three orders of 
magnitude smaller than those relative to metallic alloys. 
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These observations strongly support the scenario proposed by Bouchaud et al. [S]: 
The coalescence of the cavities is likely responsible for the observed roughness exponent 
of the C = 0.8 regime observed at larger length scales (from to 100 nm). The order 
of magnitude of is indeed in good agreement with our present observations, where 
cavities at coalescence are a few tens of nanometers wide. On the other side, the C, = 0.5 
regime observed for smaller length scales is likely related to an irregular start-stop motion 
characteristic of the non equilibrium dynamic critical-point at which the crack starts 
to advance. This may generate diffusely damped crack front corrugations naturally 
leading to a steady state roughening of fracture surfaces. This effect is predicted to 
occur even though the macroscopic velocity of the crack front is small provided that 
the instantaneous velocity during an 'avalanche' in a strongly heterogeneous medium 
may reach much higher values. This may have been at the origin of the important 
fluctuations in crack velocities reported in the present paper. 

Here, let us note that AFM observations are performed on the sample surface, 
where the mechanical state is different from that of the bulk. Hence, the observed sizes 
and growth rates of cavities at the surface may well differ from those in the bulk. New 
experiments using the FRASTA method in three dimensions applied to the post-mortem 
study of the fracture surfaces are currently being performed, in order to have access to 
the three dimensional structure of bulk damage and its evolution. Through this new set 
of experiments, one should be able to correlate also the damage structure to the fracture 
surface morphology. 

The topographic study of post-mortem surfaces will also be a good experimental 
test of the validity of Bouchaud et al. jS]. This scenario suggests that small length 
scale exponent ( = 0.5 is due to the existence of diffusively broadened crack front 
waves in localised depinning events. Other important predictions are that roughness 
amplitude should be strongly anisotropic and this effect should be very pronounced at 
low crack growth velocities. That model also predicts that the physics of the formation 
and coalescence of the cavities is responsible for the observed roughness exponent of 
( = 0.8. Such a scenario was already proved to be realistic thanks to observations of 
the roughness of growing cavities before coalescence in an aluminium alloy [7]. It will 
be of course crucial -even if much more experimentally difficult due to much lower C,c 
value- to check these hypotheses in the case of nano-ductile materials. 

Moreover, the structure of damage, which influences macroscopic mechanical 
properties such as fracture toughness and lifetime, should then be linked to the glass 
composition and nanostructure [2HI IHH tHH] • Complementary analyses addressing the 
question of the chemical bonds on the fracture surface will also be performed. 

7. Conclusion 

Quasi-static fracture in glassy materials has been studied in real time at sub-micrometric 
length scales in a wide velocity range. The crack velocities are measured as a function of 
the stress intensity factor in both amorphous materials and a partially devitrified one. 
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Important fluctuations are reported for tlie lowest velocities. They are conjectured to 
be related to the presence - inside of the glassy phase - of nano-scale damage cavities 
observed in real time ahead of the crack tip. The influence of these damage cavities 
on mechanical properties in the non purely elastic zone in the vicinity of the crack tip 
have been presented. These results confirm the scenario proposed by Bouchaud et al. |H] 
to explain the origin of the two self-affine regimes observed on fracture surfaces. The 
implications of such a ductile fracture mode in glass on the morphology of the fracture 
surfaces have been discussed. However it should be noted that our investigations are 
performed on the sample surface while the fracture surface morphology is related to the 
fracture of bulk. Consequently, it would be interesting to investigate the 3D distribution 
of damage cavities. Work in this direction is currently in progress. 

Finally, the fact that glass, at temperatures far below the glass transition 
temperature T^, joins the class of damageable materials should have important 
consequences for its mechanical properties. In applications, the design of structures 
using glass might be modified to take this behaviour into account, especially for slow 
crack propagation processes. 

The similarity between the damage modes of materials as different as glass and 
metallic alloys is an important clue to understand the origin of puzzling universal 
behaviors, hence shedding new light into the basic physical mechanisms of fracture. 
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